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Abstract 

An analytical method and a stored channel model 
were used to evaluate error performance of uncoded 
QPSK and MPSK trellis coded modulation (TCM) 
over shadowed satellite mobile channels in the pres- 
ence of phase jitter for various carrier recovery tech- 
niques. 

1 Introduction 

In coherent digital transmission systems accurate 
carrier phase and symbol timing information are 
essential for successful demodulation. The recon- 
struction of coherent reference can be accomplished 
in many ways. In general, carrier and clock are 
restored by pilot tone and self-synchronizing tech- 
niques. Self-synchronizing methods could be quite 
roughly divided into decision-directed algorithms 
and feedworward nonlinear phase estimators. This 
paper compares three representatives of the above 
classes and analyzes their relative merits on satel- 
lite mobile channels with shadowing. The impact 
of phase jitter on coherently demodulated QPSK 
and M-PSK trellis coded modulation (TCM) with 
feedforward (FF), decision directed (DD) and TTIB 
techniques for carrier phase recovery on a is consid- 
ered. 

A stored experimental channel [2] was used to 
simulate the performance of a satellite mobile com- 
munication system. The experimental data for the 
stored channel were obtained in land mobile satel- 
lite propagation measurements at L band in Syd- 


ney suburban areas via the Japanese ETS-V and 
INMARSAT-POR geostationary satellites. 

2 System Model 

2.1 System Description 

In the systems under consideration a voice codec 
generates data digits at a rate of 1/T b . If the QPSK 
modulation is used the 2-bit input sequences mod- 
ulate a QPSK waveform at a rate l/2T b symbols/s. 
In the 8-PSK TCM system the binary sequence is 
encoded in a rate 2/3 convolutional encoder. The 
3-bit symbols at the output of the encoder are then 
block interleaved and used to modulate an 8-PSK 
waveform at a rate l/2T b symbols/s according to 
the set partitioning rules [6]. The interleavear size 
is chosen on the order of the worst case fade depth 
in the simulations and assumed infinite in the anal- 
ysis. We consider three different methods for carrier 
phase recovery, FF, DD and TTIB techniques. 

2.2 Model of Shadowed Channel 

The signal in a shadowed channel is obtained as 
the sum of a lognormally distributed random phasor 
ye^° and a Rayleigh phasor : 

ae^ = ye^° + w^ y, w > 0 

where the phases 4>o and <j) are uniformly distributed 
between 0 and 2 x. 
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The conditional probability of a given y is Rician: 




where 2a represents the average scattered power 
due to multipath propagation. 

The total probability is given by 

fOO 

p(a) = / p(a | y)p(y)dy 
Jo 

It is assumed that p(y) is lognormally distributed: 


p(y) = v2tt<7i y 

with the mean 


( ( ln y-v 


0 > y > oo 
otherwise 


Py = E {y} = exp (/i + a\f 2) 


On the other hand the received sample phase can 
be expressed as 

= Oj + €\ 

where is the phase jitter and t\ is the phase error 
caused by Gaussian noise. 

In the FF algorithm for each sample we multiply 
the phase <f > j by M and perform a nonlinear trans- 
formation F on the amplitude p\ [8]. As the result 
we obtain the signal 

r’li + j r Qi = F(pi)expfjM<Ai) 

The phase of the resulting signal is 

M</>i — Ma\ + M0j = 2ki7r + M#j 

The phase error 6 is then estimated over 2N + 1 
symbols as 


and the variance 


o-y = E {(y - My) 2 } = exp(2 /x + <?l) (exp(<r?) - l) 


3 Phase Estimation in FF Tech- 
nique 

Let us assume that vj represents an MPSK symbol 
sample transmitted at time i. The phase of the ith 
transmitted sample can be written as 

a-k — 

1 1 M 

where ki is an integer (kj 6 {0,1,... ,M — 1}. 

The corresponding received sample at the input 
of the coherent demodulator is: 

rj = aiVi + nj 

where a] is a real random variable equal to the en- 
velope of the channel attenuation normalized to the 
transmitted signal, nj is a sample of a zero mean 
complex additive white Gaussian noise (AWGN). 
The recevied sample can be expressed as 

r, = r D + r Qi = piexp(j<ft) 

/v — . / r ?. -L . 


P\ = V r U + r Qi 
4>\ = tan -1 


»= (1) 
M Eil-N PiCosM(e; + < pi ) 

0 represents the phase by which the corresponding 
symbol is rotated. 

The above operations cause an M-fold ambiguity 
in the phase estimate. It can be resolved by differ- 
ential decoding for M-PSK modulation and 2 7r /A/ 
phase rotationally invariant M-PSK modulation. If 
non-rotationally invariant codes are used the phase 
ambiguity can be removed by monitoring Viterbi 
decoder metrics. 

4 Phase Estimation in DD Tech- 
nique 

The DD is a recursive algorithm which is used for 
both phase and symbol timing estimation [4]. The 
algorithm operates under the assumption that the 
incoming symbols are known to the receiver. This 
condition is met for small error probabilities or al- 
ternately if some known pattern sequence is trans- 
mitted in a training period. 

If the current phase estimate is 6\ the next symbol 
phase rotation is determined as 

0i+i = <?i + rim{rirf} (2) 

where T is a small number determined to achieve a 
trade-off between the phase error and the conver- 
gence rate and fj is the ith sample estimate. 


506 


International Mobile Satellite Conference, Ottawa, 1990 


The symbol sampling time is determined by us- 
ing a second order digital filter. The sampling time 
sequence {r,} is determined as 

n+i = n + ARe{rir*_j - r ; _irf } (3) 

where A is a small number. 

The phase error can be approximated as 

0 i+1 ~ r(ei + w) (4) 

where is the phase noise induced by Gaussian 
noise and <p\ is the phase jitter. 

In block structered data transmission the perfor- 
mance of the algorithm can be improved by sending 
a preamble at the beginning of each block [5]. 

5 Phase Estimation in TTIB 
Technique 

The filtering in the TTIB algorithm is achieved 
by using two bandpass filters with a combined fre- 
quency response equivalent to a square root raised 
cosine transfer function [7]. Each bandpass filter is 
used to isolate one of two subbands with bandwidths 
equal to one half of the symbol transmission rate. 
The subband signals are further frequency trans- 
lated to enable the insertion of the pilot tone. The 
bandstop margin is chosen large enough to accom- 
modate the worst case Doppler shift and meet the 
acquisition time requirements of the receiver. The 
reference tone is positioned centrally within the sig- 
nal bandwidth to assure an accurate estimation of 
the signal degradation. 

In the TTIB demodulator the received pilot tone 
is used to remove phase error. This process is ac- 
complished by computing the pilot tone phase and 
multiplying each symbol by the negative pilot tone 
phase. It is essential to choose an appropriate pi- 
lot tone level and filter bandwidths . As the pilot 
tone level increases more power is waisted on the 
tone transmission. On the other hand if the pilot 
tone power is reduced, it will degrade its signal- to- 
noise ratio and phase estimate. Similarly, the pilot 
tone bandwidth has an optimum. In order to reduce 
the impact of Gaussian noise it is desirable to have 
a small pilot tone bandwidth. However, the band- 
width is determined by maximum Doppler shifts in 
the channel to enable pilot tone extraction. 

The insertion of pilot tone, however, will result in 
a non-constant envelope signal and degrade the per- 
formance of this algorithm on a nonlinear channel. 


The BER degradation on a satellite channel with a 
TWT nonlinear amplifier is of the order of 0.5 dB 
[!]• 

6 Probability of Error for QPSK 

The conditional bit error probability for QPSK 
modulation given the fading attenuation due to 
shadowing a and the phase jitter 9 is given by 

P„ = Q (2£2») 

where a is the variance of the Gaussian noise in each 
signal space coordinate. 

The total bit error probability is given by 

Pb = lo I ^ p( a ) p w dad ^ 

where p(a) and p(0) are the probability distributions 
of fading attenuation and phase jitter, respectively. 

7 Error probability for TCM 
Schemes 

Let C denote the set of all possible coded signal se- 
quences {zn}, where zn is a coded symbol sequence 
of length N : 

zn = (zi,z 2 ,. . . ,zj,. . . ,zn) 

Let us assume that vn: 

VN = (vi,V 2 ,...,Vi,...,V N ) 

an element of C, is the transmitted sequence. Then 
Vj represents the MPSK symbol transmitted at time 
i. 

Vi = v D + jVQi 

where vn and vqi are I and Q channel components, 
respectively. 

For codes with large minimum number of symbols 
on error paths or for any code at small Eb/No the 
pairwise error probability can be upperbounded by 

[9] 

P(vn — » ’ Zn) = Pd,D,T 

^ 1 2<7d - r ( 1 m x \ ^ 

-2 V /4^T^ P \ 2 4 aW + ol) ^ 

where 

m x = d 2 m a m fil + Tm,tn m 
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T = E Ti H( Z Ii V Qi ~ Z Qi y Ii) 

i=l i=l 

and 

= D(m 9l 2 - rn^mj!) + T 2 (m fl22 - m\m 2 d2 ) 

N 

D = E d f 

i=l 

N 

t 2 = E T i 

i=l 

m^i = E[cos0] 
m ^]2 = Ejcos 2 #] 
m^ 2 = E[sin0] 
m ^ 2 2 — E[sin 2 0] 

The bit error probability can be upperbounded by 
Pb - kE^Ad.t (6) 

K d,D 

where ba ? D,T is the total number of erroneous infor- 
mation bits on all paths characterized by the Eu- 
clidean distance d, and the parameters D and T. 

8 Results 

The bit error rate performance with the QPSK mod- 
ulation and FF algorithm relative to the demodula- 
tion without carrier recovery on a Gaussian channel 
with variable degrees of slowly varying phase jitter is 
given in Fig. 1 The algorithm removes almost com- 
pletely the effect of phase jitter up to phase jitter 
values of 30°. The corresponding results for the 16- 
state Ungerboeck TCM code are illustrated in Fig. 
2. The algorithm operates succesfully up to phase 
jitter values of 15°. The impact of slowly varying 
phase jitter on the DD algorithm on a Gaussian 
channel is illustrated in Fig. 3. The algorithm is 
able to recover the signal phase up to 40° phase jit- 
ter. The BER for larger phase jitter exhibits high 
variations due to unreliable estimates of the signal 
phase. For a Gaussian channel with relatively low 
BER the algorithm works reliably and there is no 
need for transmiision of preambles. 

The performance of the TTIB technique depends 
on the pilot tone bandwidth. The pilot tone band- 
width should be chosen to match worst case Doppler 
shifts. For the ratio of the signal to the pilot tone 


bandwidth (3 < 50 the degradation is high. The sen- 
sitivity of the TTIB technique to phase jitter for 
QPSK and TCM is shown in Fig. 4. The TTIB al- 
gorithm is can remove slowly varying phase errors, 
but cannot eliminate the effect of Gaussian noise 
due to finite pilot tone bandwidth. The loss for the 
QPSK is about ldB for / 3 = 50. For all TTIB tests 
the power of the pilot tone level was 10 dB below 
the total output power. 

The comparative results of the three algorithms 
for the QPSK scheme on a satellite mobile shad- 
owed channel are presented in Figure 5. The per- 
formance is evaluated by the average block error 
rate (BLER) which is a suitable parameter for digi- 
tized voice transmisssion, (the block length corre- 
sponds to 30ms frames). The DD algorithm re- 
quires accurate esimates of the received signal and 
the error performance is degraded by heavy shad- 
owing. The additional information from a known 
received sequences can improve the reliability of sig- 
nal estimation. Since the voice codec generally pro- 
vides preambles at the beginning of each transmit- 
ted block of data it is possible to include this in- 
formation in the DD algorithm. It should be noted 
that most degradation in the error rate comes from 
the amplitude variations rather than phase errors, 
since phase errors in the channel were generally less 
than 10°. Therefore, carrier recovery methods could 
not achieve significant performance improvement. 
Clearly, the TTIB technique shows the superior per- 
formance. However, an additional loss of about 2dB 
should be added due to the pilot tone transmission, 
nonlinearity and Doppler shift degradations. 

The corresponding results for the 16-PSK TCM 
code is shown in Fig. 6. 

9 Conclusions 

An investigation of carrier recovery techniques on 
satellite mobile channels with shadowing has been 
performed. 

The TTIB algorithm achieves the excellent phase 
estimates for low pilot tone filter bandwidth, ac- 
commodate the worst case Doppler shifts. 

The DD algorithm with preamble transmission 
offers good performance, while the implementation 
is simple and the algorithm with the second order 
loop provides symbol timing and frequency tracking 
which is critical on mobile channels. 
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Figure 1: BER performance of QPSK with phase 
jitter on an AWGN channel: 

Dash: without FFD; 1 to 7: 0, 5, 10, 15, 20, 30 and 
40 degrees phase jitter 

Solid: with FFD; 1: 5 to 20 degrees; 2: 30 degrees; 
3: 40 degrees phase jitter 



Eb/ No ( dB) 

Figure 2: BER performance of TCM with phase 
jitter on an AWGN channel: 

Dash: without FFD, Solid: with FFD ; 1 to 5: 0, 5, 
10, 15 and 20 degrees phase jitter 
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Eb/ No ( d B) 


Figure 3: BER performance of the QPSK on an 
AWGN channel with phase jitter and DD scheme; 
Solid: Ideal QPSK; Dash: 1: phase jitter=10° 2: 
phase jitter=20° 3: phase jitter=30° 


Eb/ No ( d B) 

Figure 5: BLER performance of the QPSK modu- 
lation with the various carrier recovery schemes for 
experimental data; Dash: No carrier recovery; Solid: 
1: TTIB 2: FF 3: DD 4: DD with preamble 



Eb/ No ( d B) 


Figure 4: BER performance of the 16-state TCM 
code on an AWGN channel with phase jitter and 
TTIB scheme; Jr = 50 

Dash: No phase recovery 1: phase jitter=0° Ideal 
TCM; 2: phase jitter=30° 3: phase jitter=50° Solid: 
TTIB 1: phase jitter=30° 2: phase jitter=50° 



Figure 6: BLER performance of the 16state Unger- 
boeck TCM code carrier recovery for experimental 
data; Dash: Uncoded QPSK with ideal phase recov- 
ery; Solid: 1: TTIB 2: FF 
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